Research in Astronomy and Astrophysics, 24:045008 (12pp), 2024 April 


© 2024. National Astronomical Observatories, CAS and IOP Publishing Ltd. Printed in China and the U.K. 


https: / /doi.org/10.1088/1674-4527 /ad2ee0 


CrossMark 


Characteristics of Emission Non-nulling Pulsars Through Simulation 


Rai Yuen 


1,2,3,4 


: Xinjiang Astronomical Observatory, Chinese Academy of Sciences, Urumqi 830011, China; ryuen@xao.ac.cn 
Xinjiang Key Laboratory of Radio Astrophysics, Chinese Academy of Sciences, Urumqi 830011, China 
” Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China 
$ SIfA, School of Physics, University of Sydney, Sydney, NSW 2006, Australia 
Received 2023 November 28; revised 2024 February 16; accepted 2024 February 28; published 2024 March 27 


Abstract 


We investigate the population and several properties of radio pulsars whose emission does not null (non-nulling) 
through simulation of a large pulsar sample. Emission from a pulsar is identified as non-nulling if (i) the emission 
does not cease across the whole pulse profile, and (ii) the emission is detectable. For (i), we adopt a model for 
switching in the plasma charge density, and emission persists if the charge density is non-zero. For (ii), we assume 
that detectable emission originates from source points where it is emitted tangentially to the magnetic field-line and 
parallel to the line-of-sight. We find that pulsars exhibiting non-nulling emission possess obliquity angles with an 
average of 42°5, and almost half the samples maintain a duty cycle between 0.05 and 0.2. Furthermore, the pulsar 
population is not fixed but dependent on the obliquity angle, with the population peaking at 20°. In addition, three 
evolutionary phases are identified in the pulsar population as the obliquity angle evolves, with the majority of 
samples having an obliquity angle between 20° and 65°. Our results also suggest that emission from a pulsar may 
evolve between nulling and non-nulling during its lifetime. 
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1. Introduction 


The population and properties specific to radio pulsars that 
exhibit only non-nulling emission are unclear. Yet, pulsars with 
nulling behavior are well distinguished as a distinct type of 
pulsar (Wang et al. 2007; Konar & Deka 2019; Sheikh & 
MacDonald 2021). Identified shortly after the first pulsar was 
discovered (Backer 1970), nulling is a relatively common 
phenomenon that manifests as cessation in the whole pulse 
emission over a time that lasts from several rotation periods to 
minutes (Backer 1970; Rankin 1986; Wang et al. 2007; Young 
et al. 2015). So far, more than 200 pulsars have been reported 
with such behavior (Sheikh & MacDonald 2021). It is 
relatively straightforward to show that a radio pulsar is a 
nulling pulsar. All that is required are clear detections of 
several instances when there is no detectable emission. On the 
contrary, it is difficult to show that a pulsar is not a nulling 
pulsar (Rankin 1986), referred to as an emission non-nulling 
pulsar in this paper (or simply a non-nulling pulsar where there 
is no confusion). Observationally, there is no telling whether or 
not the thus-far continuous emission from a pulsar would be 
interrupted in subsequent pulses. One difficulty lies with the 
similarities in the emission properties (when not nulling) and 
the rotation characteristics between the two types of pulsars. 
Therefore, it is challenging to predict whether a pulsar will 
exhibit non-nulling emission based only on its rotation or 
emission characteristics. Nevertheless, there were many 
investigations into the properties of nulling, which also offered 


hints on certain properties for pulsars that should not exhibit 
nulling. For example, the suggestion that null pulses appear to 
increase with the pulsar’s age also implies that pulsars with 
large age are less likely to exhibit non-nulling. Ritchings 
(1976) concluded that nulling is related to the rotation period of 
a pulsar, such that pulsars with longer periods should null more 
often (Biggs 1992; Wang et al. 2007). There were also studies 
suggesting that half of the known pulsars should demonstrate 
nulling, which also imply that half of the known pulsars may be 
non-nulling. There were predictions that the frequency of null 
is related to the morphology of the integrated profiles 
(Biggs 1992). By separating profiles into different classes, 
Rankin (1986) showed that pulsars with pulse profiles having a 
single component should null less frequently (Weisberg et al. 
1986) than pulsars possessing profiles with multiple compo- 
nents. Biggs (1992) further suggested that a pulsar would null 
more when its obliquity angle, a, evaluated between the 
rotation and magnetic axes, is small. However, a statistical 
study of nulling pulsars concluded that nulling shows less 
correlation with the obliquity angle or the rotation period but 
more with the profile width (Li & Wang 1995). This indicates 
that nulling, and non-nulling, may also be related to the 
emission geometry. 

The purpose of this paper is to investigate several 
characteristics of emission non-nulling pulsars through studies 
of the distribution of a, 8 and different profile widths, or duty 
cycles, 6, from simulation of a large pulsar sample. Here, 
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G=C-—a is the impact parameter, and Ç is the viewing angle 
between the rotation axis and the line of sight (Radhakrishnan 
& Cooke 1969). There are two criteria for confirming a radio 
pulsar as an emission non-nulling pulsar, namely (i) its 
emission does not cease or cessation does not occur across 
the whole pulse profile, and (ii) its emission must be detectable. 
Point (i) involves the knowledge for changes in the pulsar radio 
pulses between emission (non-nulling) and nulling. An 
important source of such information comes from the 
consideration of the quasi-periodic emission switching in 
intermittent pulsars between pulse “on,” when emission is 
detectable, and “off,” in which the emission ceases (Kramer 
et al. 2006; Camilo et al. 2012; Lorimer et al. 2012; Lyne et al. 
2017). The mechanism of the phenomenon can be interpreted 
as changes in the emitting plasma density between two states of 
ample and vacuum (Kramer et al. 2006). This implies 
correlation between radio emission and the plasma density, 
and pulse cessation arises when the plasma density is zero. The 
correlation is consistent with the common view that pulsar 
radio emission is produced from ultra-relativistic plasma out- 
streaming along the open field-lines in two-stream instability, 
which leads to the prediction that the pulse intensity is 
proportional to the plasma density (Manchester & Taylor 1977; 
Cordes 1979; Lyubarskii 1996). Although there is still no 
commonly accepted model for changes in the plasma density, a 
model for multiple quasi-stable emission states, each with 
different plasma charge densities, has recently been proposed 
for obliquely rotating pulsar magnetospheres (Melrose & 
Yuen 2014). The model is based on synthesis of two limiting 
conditions for plasma charge density. One limit corresponds to 
a corotating magnetosphere, in which the plasma charge 
density maintains the Goldreich & Julian (1969) value, poy. 
Another limit refers to the vacuum-dipole model, in which the 
origin of the inductive electric field, Eina, is attributed to a 
rotating magnetic dipole. The need to screen the parallel 
component of Eina requires the presence of a charge density, 
denoted by p,,, that is arbitrarily small and insufficient to 
generate the potential field required by corotation. A class of 
synthesized models for emission states characterized by the 
parameter y is then constructed through consideration of the 
plasma charge density in each emission state as a combination 
of y and (1—y) fractions of p,;, and poy, respectively. 
Consequently, a change in the value of y would result in a 
change of the plasma charge density, which corresponds to a 
change in the emission state. It follows that the emission will 
persist in any emission states in which the plasma charge 
density is non-zero. In order for the changes in pulse emission 
to be detectable, the switching in the emission state must occur 
at the locations that are visible to an observer. Here, we adopt 
an emission geometry in which radiation is emitted from source 
points on the field-lines whose footpoints are located in the 
open-field region. Visible emission comes only from the source 
points (visible points) where the emission is directed 


tangentially to the dipolar magnetic field-lines (Cordes 1978; 
Hibschman & Arons 2001; Kijak & Gil 2003) and parallel to 
the line-of-sight direction. The location of a visible point in the 
magnetosphere for a pulsar phase, p, can be determined 
uniquely given ¢ and a. As the pulsar rotates, the visible point 
moves tracing a closed path for every pulsar rotation, referred 
to here as the trajectory of the visible point. Changes in the 
pulse emission can be detected by an observer provided that the 
changes take place at the locations defined by the trajectory of 
the visible point. 

The observations of more and more pulsars with detectable 
switching in the emission properties (Smits et al. 2005; Kramer 
et al. 2006; Rankin et al. 2006; Lyne et al. 2010, 2017; Camilo 
et al. 2012; Lorimer et al. 2012; Keith et al. 2013; Marshall 
et al. 2015; McSweeney et al. 2017) demonstrate that it is 
common for pulsar magnetospheres to undergo sudden changes 
of some sort. An essential piece of information from these 
observations is that the timescale for emission switching in 
different pulsars is different. It follows that the switching 
timescale in some pulsars may be too long to be captured in 
regular observing sessions. For our investigation, we assume 
that emission switching occurs in all our samples. From study 
of the polarization position angles based on a large sample of 
pulsars (Rankin 1993a; Tauris & Manchester 1998), the 
distribution of œ was demonstrated to be situated within the 
range of a < 90°, with the distribution occupying mostly small 
values (Rankin 1990). The result was reinforced by different 
simulations (Zhang et al. 2003; Kolonko et al. 2004). Hence for 
our samples in this paper, we select a values from the range 
between 1° and 90°. In addition, the width of the profile must 
be known to determine the emission across the profile. We 
consider duty cycles from 0.05 to 1 obtained from different 
pairs of Ç and œ under an assumed maximum emission height. 
As for the last case, pulsar radio emission was shown to 
originate from a height at about 0.17, (Johnston & Weis- 
berg 2006), with r, = c/w, being the light-cylinder radius and 
w, the spin frequency of the star. In our simulation, a maximum 
height of 0.27, is assumed. We also focus on emission from the 
main pulse without referring to the interpulse. Some pulsars 
were observed to demonstrate correlated changes between 
different emission features, such as subpulse drifting, nulling, 
quasi-periodic modulation, micro-structure and polarization 
properties (Taylor et al. 1975; Bartel et al. 1982; Rankin 1986; 
van Leeuwen et al. 2002; Janssen & van Leeuwen 2004; 
Redman et al. 2005), suggesting that variations in the radio 
emission may involve a more complex change in the emission 
process. Despite many investigations, such correlation between 
different emission phenomena is still largely uncertain and its 
relationship to the pulsar radio generation process remains 
unclear. There are also unconventional models for the pulsar 
emission process which involve unusual magnetic field 
structures (Gao et al. 2016, 2017) and a pulsar model with 
strange quark matter (Xu et al. 1999, 2001). Before including 
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such complicated cases, it is useful to explore the implications 
of the traditional definition in an idealized model. 

We outline the model for non-nulling emission in Section 2. 
The simulation procedures and the results are presented in 
Section 3. The analysis for the correlations between emission 
non-nulling pulsars and the pulsar parameters is shown in 
Section 4, and we discuss and conclude the paper in Section 5. 
The electromagnetic fields, the emission geometry and the 
transformation matrices referred to in the paper are outlined in 
Appendices A-C, respectively. 


2. Model for Non-nulling Emission 


Verifying the emission of a pulsar as non-nulling involves 
ensuring that the visible emission does not cease across the 
whole pulse profile in every rotation. In this section, we 
summarize the model presented by Melrose & Yuen (2014) 
with an emphasis on changes in plasma density in relation to 
the condition for non-nulling emission. 


2.1. Detectable Changes in the Plasma Density 


The interpretation of emission changes in intermittent pulsars 
as due to changes between two states of plasma density 
(Kramer et al. 2006) implies that different emission states, each 
with unique plasma density, are allowed in a pulsar magneto- 
sphere. In the model for pulsar magnetospheres with multiple 
emission states (Melrose & Yuen 2014), the different emission 
states can each be identified by a particular value in the 
parameter y between two limiting conditions, designated by 0 
and 1, respectively. The two conditions can be defined through 
the electric field in an obliquely rotating pulsar magnetosphere, 
which takes the form 


E=( —ybb)- Ena + (1 — y)Epot, (1) 


where b is the unit vector along the dipolar field-lines. In 
Equation (1), Eina signifies the electric field induced by an 
obliquely rotating magnetic dipole, with the form given by 
Equation (A2), and the electric field associated with the 
corotation charge density is given by E,o.=—V cor. The 
limiting condition for y= 0 then corresponds to the corotation 
electric field, E = Eeor, given by Equation (A5) as specified in 
the corotation model (Goldreich & Julian 1969), and the 
associated corotation charge density is signified by 


Pol = —EqV $ Eor- (2) 


The other limit, represented by y= 1, is associated with the 
“minimal” model, which is signified by the screening of the 
parallel component of the inductive electric field, Eina, while 
the perpendicular component, Einaı, has the same value as in 
the vacuum-dipole model. The screening requires a charge 
density, Pin, to generate an electric field such that 


Emin = —Ejna|. The p,,,, takes the form given by 
Pmin = — £0 V i (b Einaj). (3) 


The necessity to screen the electric field parallel to the magnetic 
field-lines in every emission state for an oblique rotator 
demands a charge density, which has the form given by 
Melrose & Yuen (2014) 


Psn = YPmin T (a — Y) Por (4) 


The screening is effective if the condition is met for the 
creation of the secondary particle pairs that result in a 
multiplicity, À, larger than unity (Beskin et al. 1993; Melrose 
& Yuen 2012). The last case can be written as 
A= e(n, +n_)/pcy>>1 (Gurevich & Istomin 2007), with 
Poy = e(n, — n_) and n4 being the number density of positrons 
and electrons. In our model, the plasma density is signified by 
Ap,,/€. The correlation between the pulse intensity and the 
plasma density (Manchester & Taylor 1977; Cordes 1979; 
Lyubarskii 1996) implies that the former is proportional to 
APs, /€. Assuming the same A across different emission states 
would mean that the pulse intensity is proportional to p,,. It 
follows that a variation in the value of y corresponds to a 
change in p,,, which results in a change of the pulse intensity. 

The value of y, and the associated value of p,,, at given 
coordinates in a magnetosphere can be determined uniquely in 
an emission geometry defined by a pair of ¢ and a. We adopt 
an emission geometry (see Appendix B) in which the magnetic 
field-lines are of a dipolar structure, and emission arises from 
source points on the field-lines whose footpoints are located 
within the open-field region (Cordes 1978; Kijak & Gil 2003). 
The condition for visible emission then requires that the 
emission be emitted tangentially to the local field-lines (J. A. 
Hibschman & J. Arons 2001) and in parallel to the line-of- 
sight. From such geometry, an explicit solution can be obtained 
for the source point of visible emission in terms of the polar 
and azimuthal angles as a function of 7. In the magnetic frame, 
the angles are represented by 6,y and py and their expressions 
are given by Equation (B1). Transforming to the observer’s 
frame using Equation (C5), the visible point can be expressed 
in terms of the angles 0y and ¢y relative to the rotation axis. 
For an oblique rotator, the visible point moves forming a closed 
path after one pulsar rotation, which is called the trajectory of 
the visible point. Figure 1 shows three examples for the 
trajectory of the visible point in the magnetic frame. Note that 
the size and shape of a trajectory vary depending on the values 
of ¢ and a, and they are generally not circular and may or may 
not include the magnetic pole. 

The assumption that emission originates from the open-field 
region implies the dependence of the visible point on height. 
Given a pair of Ç and a, there exists a minimum height, ry, at 
which the trajectory of the visible point is tangent to the locus 
of the last closed field-lines (see Appendix B). The height, ry, 
changes as a function of w as indicated by Equation (B2). To 
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Figure 1. Plot illustrating three trajectories of the visible point using a = 25° 
and ¢ = 10° (green), 30° (brown) and 60° (blue). The plot is constructed in the 
magnetic frame such that the magnetic pole is located at the origin. Both blue 
and brown trajectories enclose the magnetic pole. The boundary of an open- 
field region at 0.27, is drawn in black dotted curve with its center located at the 
origin. Also shown is the non-nulling (gray) region for the same a. While all 
the trajectories cut the open-field region differently, only the brown and green 
trajectories traverse the non-nulling region, with the latter trajectory lying 
entirely inside both the open-field and non-nulling regions. 


illustrate that, Figure | also shows an open-field region with a 
boundary (black dotted closed curve) at 0.2r,. All trajectories 
cut the open-field region but at different pulsar phases. For the 
blue trajectory, it enters the region at A, corresponding to 
pı = —30°, and exits it at B, where 7 = 30°. The rest of the 
trajectory is located at heights that are above 0.2r,. Similarly, 
only the part of the brown trajectory between C and F lies 
inside the open-field region, with %1 2 = 94°, respectively, 
between which emission is visible. This is different for the 
trajectory in green, which lies entirely within the open-field 
region, meaning that the variation of height along the trajectory 
is always below 0.2r,, and emission is potentially visible 
throughout the whole pulsar rotation. The range of pulsar phase 
for detectable emission is bounded by p; and 72, and hence the 
profile width is given by Aw = Y2 — Y1, or equivalently in duty 
cycle by 6= Ay /2z. For the blue trajectory, Ay = 60°, or 
60.17, whereas it is Ay = 188°, or 60.52, for the brown 
trajectory. The green trajectory, being completely enclosed in 
the open-field, has a Ay = 360°, which gives ô= 1. 


2.2. Geometry for Non-nulling Emission 


We define non-nulling emission as visible emission that does 
not cease across the whole pulse profile in every pulsar 
rotation. For a pair of ¢ and a, the angular locations of the 
visible point can be uniquely defined by the trajectory of the 
visible point. Then, Equation (4) is used to solve for the y value 
that satisfies p,, = 0 for each Y = Yg between Yı and y» that 


lies inside the open-field region. A pulsar sample is regarded as 
an emission non-nulling pulsar provided that no solution exists 
for p,, = 0 with O < y < 1 at one or more wx. This implies that 
the open-field region and the non-nulling region must contain 
overlaps and the overlapping regions must be cut by the 
trajectory of the visible point. Such relationship between the 
emission geometry and non-nulling emission is illustrated in 
Figure 1. The figure demonstrates a region (in gray) in which 
Psy > O for O<y< 1. It is clear from the figure that the 
location for the non-nulling region is different from that for the 
open-field region. This implies that a trajectory that cuts the 
latter does not necessarily also cut the former. An example is 
the blue trajectory, with which it does not cut any part of the 
non-nulling region even though it lies partially in the open-field 
region. It follows that the emission along the blue trajectory 
does not satisfy our definition for non-nulling emission because 
it lies entirely in the region where p,, = 0 exists for 0< y< 1 
(cessation across the whole profile). The brown trajectory 
traverses the open-field region from C to F, and cuts the nulling 
region only from D to E at Yı 2 = F 65°5. This means that the 
parts of the trajectory through CD and EF will always be 
emitting (cessation not occurring across the whole profile), and 
hence the pulsar would be regarded as an emission non-nulling 
pulsar. With ô= 1, the green trajectory is immersed completely 
in the open-field region, and visible emission occurs along the 
entire trajectory which also lies within the non-nulling region 
giving non-nulling emission detectable across the whole 
profile. 


3. Simulation Procedures and Results 


Our investigation of emission non-nulling pulsars and their 
correlations with the duty cycle (6), the impact parameter (3) and 
the obliquity angle (œ) is based on simulation of a large sample 
of pulsars. Below, we summarize the simulation procedures in 
relation to the criteria for identifying non-nulling emission. 

In three-dimensional coordinates, the rotation axis of a pulsar 
and the line of sight of an observer are both random and 
uniform resulting in an isotropic distribution. We adopt a 
random distribution of the obliquity angle described by a 
cosine function (Smits et al. 2009). Furthermore, an emission 
from a pulsar is assumed to be beamed toward an observer, 
meaning that the values of @ are dependent on the distribution 
of a. We assume the likelihood of the magnetic axis occurring 
in the direction of the line of sight being proportional to 1/47. 
For our simulation, we consider the ranges of 1° < a < 90° and 
0° < ¢< 90°. We also restrict 3 to within the value of |6| < 15° 
for consistency with observations (Lyne & Manchester 1988), 
and apply an upper limit of r < 0.2 r for the emission height to 
all the samples (Johnston & Weisberg 2006). We assume that 
all the sample pulsars are emitting in an emission state with 
0<y <1, and each is capable of emission state switching in 
association with changes in p,,. With these parameters defined, 


Research in Astronomy and Astrophysics, 24:045008 (12pp), 2024 April Yuen 


1.0 


0.8 


0.6 


Counts 


0.4 


0.2 


0.0 
0.0 0.2 0.4 0.6 0.8 1.0 


6 


Figure 2. Histogram showing the normalized distributions of emission non- 
nulling (blue) relative to nulling (gray) pulsars for the 20 values of 6 between 
0.05 and 1.0. 


the non-nulling emission of a sample with a pair of ¢ and a is 
determined at {6y, dy} for each w from —180° to 180°, in steps 
of O°1, along the trajectory of the visible point using 
Equations (B1) and (C5). The values of vw, and y) are 
determined, making sure that the heights are below 0.27, with 
the fiducial plane assumed to be at 7) = 0°. After that, violation 
of the condition for p,, = 0 with O < y < 1 is checked for each 
wg, between Y; and %2, using Equation (4). The pair Ç and a is 
labeled as non-nulling only if such violation is identified. We 
obtain a total of 357,395 samples for non-nulling and nulling 
pulsars. 

Figure 2 shows the results of the distribution for the emission 
non-nulling pulsars for different 6, in steps of 0.05, from the 
simulation. For comparison, the corresponding distribution for 
nulling pulsars is also shown. We find that emission non-nulling 
pulsars exist for all 6 from 0 to 1. Specifically, the proportion of 
non-nulling pulsars increases relative to that of the nulling 
pulsars as 6 increases. The former dominates from 6 > 0.2, and it 
reaches 100% from 6=0.4. We obtain 289,441 detectable 
samples, giving 81.0% of non-nulling pulsars. Overall, about 
53% of the non-nulling pulsars have 6 < 0.2, and around 80% 
possess a duty cycle of <0.45. Comparing with the estimated 
6~0.1 for ordinary pulsars (Kolonko et al. 2004; Maciesiak 
et al. 2011), our result indicates that the distribution of emission 
non-nulling pulsars contains relatively large duty cycles. 


4. Properties of Non-nulling Pulsars 


In this section, we examine the relationships between non- 
nulling emission and each of (3, a and ô, and the implications of 
the correlations. 


4.1. Relations to the Emission Geometry 


The variation in the distribution of a for emission non- 
nulling pulsars for ten different 6 values is shown in Figure 3. 
As 6 increases, the total number of non-nulling pulsars changes 
in such a way that it rises and reaches a maximum at 6 ~ 0.2, 
followed by a decrease and arriving at the minimum at 6 = 1.0. 
At small 6, the distribution covers all the œ values in our 
simulation. Beyond that, the number drops and reaches zero at 
a= 90°. As ô increases, the distribution undergoes significant 
changes from spreading over all the a values to increasingly 
shrinking to cover only a smaller range of œ up to ~25° at 
6=1.0. Furthermore, non-nulling pulsars are not found for 
a = 30° from 6 = 0.4, whereas they always exist at a < 30° for 
all 6. In addition, the total number under the distribution 
decreases considerably. A peak showing the most common a 
value is also clearly seen in each distribution, which shifts 
toward smaller a as ô increases. The last is consistent with the 
mean value of a that reduces as 6 increases. For 6 > 0.5, the 
changes are mostly seen as a decrease in the total pulsar 
number and shrinking of the distribution, but the general shape 
of the distribution remains. This suggests a weaker dependence 
of the distribution on large 6. 

An alternative presentation for the variation in the number of 
non-nulling pulsars as a function of 6 for different œ values is 
given in Figure 4. The a values are divided into nine 
consecutive intervals. For a between [0, 10)°, the pulsar 
number increases as 6 increases from 0.05 to about 0.4, beyond 
which the number is almost constant. For a values up to about 
20°, the number of non-nulling pulsars at large 6 (>0.4) is more 
than that at small 6. This reverses as the values in an a interval 
increase, and the curves display a clear maximum that shifts 
increasingly toward smaller 6. Beyond the maximum, the 
pulsar number drops, but at a rate that is unique for different a 
intervals. The 6 range under a curve reduces gradually as the 
values in an qa interval increase. It begins with a non-zero 
pulsar number over all 6 values (brown, blue and gray) to reach 
zero at a 6 value that decreases as a increases. Note that almost 
similar changes occur for the two a intervals from 70° to 90°. It 
is apparent from Figures 3 and 4 that the number of emission 
non-nulling pulsars is dependent on both 6 and a. 

The variation in the distribution of 3 is shown in Figure 5 for 
the same values of 6 as in Figure 3. For all 6 in our simulation, 
we find that there are always non-nulling pulsars that possess 
B <0°, with the amount being dependent on 6. However, the 
number of non-nulling cases with positive ( increases as 6 
increases, but both +6 and —£ are never distributed evenly 
regardless of 6. This can also be seen from the mean value of 3 
that displays increasingly positive behavior as 6 increases in 
general, but never equals zero. Furthermore, the total number 
of pulsars is largest for 6 ~ 0.2 and decreases as 6 increases. 
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Normalized counts 


Figure 3. Histogram showing the normalized distribution of a for ten different 6. The 6 value increases from left to right and from top to bottom. The range of a under 
the distribution reduces as 6 increases. The mean value of a for each distribution, @, is given for each ô. 


4.2. Evolution of the Correlation between Non-nulling 
and a 


As displayed in Figure 3, the changes of covering all a at 
small 6 to being increasingly concentrated around smaller a as 
6 increases imply that emission non-nulling pulsars with 
different 6 show different preferences for a values. An 
alternative presentation is depicted in Figure 6, which shows 
a normalized distribution for the number of non-nulling pulsars 
as a function of a. Noticeable fluctuation is seen in the number 


of non-nulling pulsars across the a values, but an overall 
pattern can still be recognized. In general, emission non-nulling 
pulsars can be found for almost any a values in our simulation. 
However, the distribution varies in such a way that it possesses 
a certain number of non-nulling pulsars at œ= 1°, which 
increases linearly as œ increases and reaches a peak at a ~ 20°. 
After that, the distribution decreases approximately exponen- 
tially as œ increases up to about 65°, followed by a roughly 
constant distribution for the next 20°. It then drops again 
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Figure 4. Plot showing variation of the normalized pulsar number counts for different intervals of a as a function of 6. The legend on the right indicates the nine 
intervals of a with a square bracket, and a parenthesis signifies the corresponding boundary value being included or excluded, respectively. Note that the curves for the 


two largest intervals are overlapped for 6 = 0.2. 


rapidly reaching zero at a=90°. We find that 50% of the 
pulsars have a < 39°5, and 80% possess a < 66°, with the 
average for a being 42°5. This demonstrates that emission non- 
nulling pulsars tend to possess smaller œ values, which is 
consistent with observations (Tauris & Manchester 1998). Also 
shown in Figure 6 is the normalized cumulative sum, as a green 
curve, representing a function of increasing a. In general, the 
slope of the curve is positive meaning that the cumulative count 
increases as a increases. However, the slope is not constant, 
and a rapid increase in number occurs at small œ up to about 
30°, followed by a slower increase as a increases further. 

To adequately compare our results with observations would 
require separating the non-nulling pulsars from their nulling 
counterparts. As mentioned earlier, it is difficult to confirm 
from observational data alone that the emission from a 
currently emitting pulsar will never cease. This means that 
estimation of a for pulsars in any survey or sample will 
inevitably mix together the two types of pulsar. However, in 
addition to previous investigations that showed pulsar obliquity 
angles are small (Tauris & Manchester 1998), recent measure- 
ments of œ from 190 new Five-hundred-meter Aperture 
Spherical Telescope (FAST) pulsars (Wang et al. 2023) have 
also demonstrated similar results. Using the rotating vector 
model (Radhakrishnan & Cooke 1969), the a values of these 
pulsars are reported between 0° and 180°. Since the field-line 
structure is symmetric between the northern and southern 
magnetic hemispheres for a dipolar magnetic field, this gives an 
average of a ~ 37° (from the near rotation axis) and a ~ 140° 
(or ~40° measured from the far rotation axis), with || most 
likely to be <5°. Again, it is likely that the pulsar sample 


contains both non-nulling and nulling pulsars. Considering that 
nulling pulsars tend to have small obliquity angles (Ritch- 
ings 1976; Cordes & Shannon 2008), and that about 20% of the 
samples in our simulation are nulling pulsars, we remove 20% 
of the pulsars with the lowest a values from the 190 pulsars 
(Wang et al. 2023). The average a values from the remaining 
pulsars are about 44° and 133° (or 47° from the far rotation 
axis), which are consistent with our results. 

There are suggestions for correlation between pulsar age and 
the evolution of the obliquity angle. From the consideration of 
energy loss through magnetic-dipole radiation of an oblique 
rotator (Manchester & Taylor 1977; Lyubarskii & Kirk 2001), 
pulsar braking decreases as a decreases. The braking reaches 
zero at œ = 0°, where there is no electromagnetic radiation. This 
implies that evolution of a is toward smaller values as a pulsar 
ages. Based on this trend for changes in the a value, Figure 6 
may be broadly divided into three phases each representing 
different evolutionary properties in the number of emission non- 
nulling pulsars. The first phase (J) is signified by the nearly 
constant number of emission non-nulling pulsars across the 
range from about 90° — 65°, in which little correlation is seen 
with the œ values. As a decreases further, the number of non- 
nulling pulsars shows a near exponential increase that ends at 
a 20° in the second phase (ID. The last phase (Z) covers 
a < 20°, during which the correlation between œ and the number 
of pulsars with non-nulling emission is reversed such that the 
latter decreases linearly as the former decreases. We find that the 
majority of the pulsars lies in phase J giving a proportion of 
57.5% of the whole sample, and leans more toward small a. The 
next highest number of pulsars is found in phase J, which 
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Normalized counts 


Figure 5. Similar to Figure 3 but for the normalized distribution of 8 showing their distributions that cover an increasingly broader range of 8 as ô increases. The mean 


value of G is given for each ô. 


contains about 22.2% of the sample. The pulsars in phase Z 
constitute about 20.3% of all the samples. The drop in the 
number of non-nulling pulsars at large a leads to the conclusion 
that lesser younger pulsars should exhibit non-nulling emission. 
In our model, this is due to fewer and fewer pulsars with both 
large a and ô> 0.3, as shown in Figure 4. From the geometry 
outlined in Section 2, the scenario in which emission comes 
from a broad 6 while keeping the source points at lower heights 
is increasingly rare as œ increases. As pulsar radio emission is 


assumed to come from the inner magnetosphere, this results in 
the decrease in the number of non-nulling pulsars (as well as the 
overall detectable pulsars) with large a. The continual decrease 
in the pulsar number as a decreases in phase J implies that lesser 
and lesser pulsars demonstrate non-nulling emission as they age. 
The existence of variation in the number of non-nulling pulsars 
across different œ values suggests that non-nulling emission in 
the pulsars is not fixed, but it changes over the lifetime of a 
pulsar. 
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Figure 6. Plot showing variation in the normalized number of emission non-nulling pulsars as a function of a as a blue curve. The green curve represents the 
corresponding cumulative sum normalized to its maximum value. The boundaries for three phases along the blue curve are J := [65°, 90°], Z := [20°, 65°) and IH := 


[1°, 20°). 


5. Summary and Discussion 


We have identified several characteristics of emission non- 
nulling pulsars. From the distributions of obliquity angle, we 
discovered that the total number of non-nulling pulsars exhibits 
significant decrease as 6 increases up to 0.5, along with 
decreases in the range of œ toward covering smaller values. 
However, the distribution is roughly consistent for 6 > 0.5. By 
assuming a decreases as a pulsar ages, we identified three 
evolutionary phases for change in the pulsar number. Almost 
half of the pulsars favor the evolutionary phase signified by a 
between 20° and 65°, in which the pulsar number increases as 
a decreases. We also found slightly more non-nulling pulsars 
with —@ than those with +68 for large 6, but they are 
increasingly biased toward —(@ as the 6 decreases. One 
prediction of the model would be that pulsars with large 6 
and small œ are likely to exhibit non-nulling emission. An 
example relates to PSR B0826—34 and the confirmation of a 
weak emission mode in the pulsar (Esamdin et al. 2005). 
Emission from the pulsar can be detected in both the main 
pulse and interpulse, with each having a broad profile window 
of about 100°. The pulsar is estimated to have a small a of 
about 0°5 (Esamdin et al. 2005), and an equally small 8 with a 
positive value (Lyne & Manchester 1988; Rankin 1993b; 
Gupta et al. 2004). Such conditions fit well with our criteria for 
non-nulling emission, and hence the pulsar is predicted to be a 
non-nulling pulsar in our model. This is consistent with the 
latest observation (Esamdin et al. 2005) that confirmed the 
previously thought null emission (Biggs et al. 1985) was in fact 
a weak emission mode. Therefore, the pulsar is an emission 
non-nulling pulsar. For ordinary radio pulsars with estimated 


average profile width of 6~ 0.1 (Manchester & Taylor 1977; 
Maciesiak et al. 2011), our model predicts that roughly half 
would be emission non-nulling pulsars. 

It is apparent from Figure | that emission non-nulling pulsars 
and observed emitting pulsars are not the same. The former 
would be represented by the pulsar with the green trajectory, 
and the latter would possess a trajectory that cuts an 
overlapping region from the nulling and the open-field regions 
similar to the blue trajectory. In our simulation, emission state 
switching is assumed in all pulsar samples without the need to 
specify the frequency of switching (or simply the switch rate). 
This is because our definition for non-nulling emission requires 
a part or whole of the trajectory of the visible point to lie 
simultaneously within the non-nulling region and the open- 
field region. This means that the emission from a non-nulling 
pulsar will never cease regardless of emission switching and 
the switch rate. This is the case for the green trajectory. 
However, different switch rates have different consequences on 
classifying a pulsar whose visible emission comes from the 
trajectory of the visible point that lies in the nulling region. 
Here, a low switch rate would mean that the wait time for 
switching of p,, between zero and non-zero is long, and short 
for a high switch rate. It follows from the discussion in 
Section 2.1 that pulsars with different switch rates to p,, = 0 
would result in different probabilities for detecting the null 
pulses. Consider the blue trajectory as an example, which 
traverses both the open-field region and the nulling region 
(from A to B) and pulse cessation will take place along the part 
of the trajectory when switching to p,, = 0 occurs in the latter 
region. For a high switch rate, an emission state switching of 
Psn between non-zero and zero is frequent, and the pulsar is 
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likely reported as a nulling pulsar. For a low switch rate, there 
exist two conditions for p,,. Assume first the pulsar is initially 
in the state of p,, + 0. A low switch rate means that switching 
to Pn = O is infrequent and could escape the detection from 
different observations, and so the pulsar would be reported as a 
non-nulling pulsar. This has the consequence of low detection 
rate for nulling pulsars even though the results of our 
simulation for 6~0.1 suggest that non-nulling and nulling 
pulsars each constitute nearly half of the pulsar population. 
Alternatively, the pulsar may be initially in an emission state 
where p,, = 0. The wait time is long for switching to another 
emission state in which p, + 0, and detecting the emission 
from the pulsar would require multiple and long observations. 
The latter case has implications on searching for new pulsars. 
Pulsar searching involves collecting signals from a pointing in 
the sky, and the time spent on each pointing is usually limited. 
For example, the average observation time for each pointing as 
outlined in the Parkes Pulsar Survey was 157.3 s (Manchester 
et al. 1996). A newer survey, known as the FAST Galactic 
Plane Pulsar Snapshot (Han et al. 2021), which discovered 
about 200 new pulsars, was performed using an integration 
time of about 300s for each pointing. With these observing 
durations, discovering emission changes from pulsars with low 
switch rates may be difficult, implying that a certain number of 
pulsars will inevitably be missed in a survey. The above 
suggests that a true pulsar population with accurately 
distinguished pulsar types would require observations with 
long duration so that pulsars exhibiting emission state switch- 
ing between different switch rates can be adequately captured. 

Our results show that the emission status (non-nulling or 
nulling) of a pulsar may evolve over time. The assumption of 
radio pulsars possessing random values of Ç and a implies that 
the 8 values should be equally distributed between negative 
and positive. Let us assume that a pulsar initially possesses a 
negative ĝ, and it is an emission non-nulling pulsar, similar to 
that illustrated in Figure 1. As the pulsar ages, the value of a 
decreases, and the 3 becomes more positive. From Figure 1, the 
latter implies that the trajectory of the visible point would 
slowly enter the nulling region. On the other hand, a pulsar that 
begins with a positive @ will remain positive. Therefore, the 
change in the visible emission is toward nulling in both cases. 
However, if the change of a is from small to large, as suggested 
from energy loss in relation to the longitudinal current flow in 
the magnetosphere (Beskin et al. 1988), a positive © would 
become more negative as a pulsar ages. In this case, the 
trajectory of the visible point will likely leave the nulling 
region at some point in time. Here, the change in the visible 
emission is toward non-nulling. For a pulsar with the trajectory 
in the nulling region, it may or may not exhibit emission 
switching that leads to pulse cessation as described in 
Section 2.1. However, observations have found more and 
more pulsars with the properties of emission switching (Smits 
et al. 2005; Kramer et al. 2006; Lyne et al. 2010, 2017; Camilo 
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et al. 2012; Lorimer et al. 2012; Keith et al. 2013; Marshall 
et al. 2015). This strongly suggests that changes of some kind 
in the magnetospheres are common and that their numbers are 
increasing among radio pulsars. This is consistent with more 
than 200 currently known nulling pulsars (Sheikh & MacDo- 
nald 2021), which indicates that the phenomenon is relatively 
common. Regardless of the model for pulsar energy loss, the 
examination for the evolution of radio pulsars based on the 
above predictions would require high quality data to accurately 
measure the values of ¢ and a. In addition, our definition for 
non-nulling emission is based on the assumption that emission 
takes place for p,, > 0 and emission switching to p,, = 0 does 
not occur across the whole profile width. It means that pulsars 
with the value of p,, close to zero would still be classified as 
non-nulling pulsars. Detection of such presumably weak 
emission would require a powerful telescope. The latter is 
realized given the observing power of FAST (Li et al. 2018), 
and the future availability of telescopes and large arrays, such 
as the Qitai Radio Telescope (QTT) and the Square Kilometre 
Array (SKA). With the prediction of discovering more than 
2000 new pulsars (Xie et al. 2022), the 110m QTT 
(Wang 2014) is promising to provide higher quality data for 
revealing radio emission of increasing details, thus enhancing 
our investigations of the pulsar radio emission mechanism. 
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Appendix A 
The Electromagnetic Fields 


Considering an obliquely rotating magnetic dipole, js, the 
time-dependent vector potential, A, is of the form (Melrose & 
Yuen 2014) 


(Al) 
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where the position vector from the stellar center is represented 
by x, with r= |x|. The associated inductive electric field is 
given by 


Bge H(t) (A2) 


4T r? 


When expressed in spherical coordinates, Equation (A2) can be 
written in the form 


Eind, mw sin a 0 
fe | L cos( — y) 


cos 0 sin(ġ — p) 


Eina, | = (A3) 


2 
Tr 
Eind,¢ 
The Eina is proportional to sina, and so it vanishes in an 
aligned rotator. The corresponding magnetic field equation is 
given by 


An 5 4 


uP a pr 
p= tol Be H rh, 3xx u-— rù 
F rE 


+ 


EXE] (A4) 


r 

where the radiative terms are expressed as œ1/ r and xl fr, 
and the dipolar term is signified by œ1/ r°. Note that p is a 
function of œ in the observer’s frame, and Epot: Eina and B are 
all functions of a. 

The electric field vanishes in a plasma-filled magnetosphere 
with negligible particle inertia and infinite conductivity, which 
gives 


Eor = — (w, X x) X Biip» (A5) 


in the co-moving frame, where w, is spin frequency of the star. 
Equation (A5) represents the corotation electric field (Gold- 
reich & Julian 1969). For a magnetosphere in oblique rotation, 
Equation (A5) may be written as (Hones & Bergeson 1965; 
Melrose 1967) 


Egor = —V cor >. A 

Ot 
where Eina = —OA/Ot. When expressed in spherical coordi- 
nates, Ecor has only the radial and polar components and is 


perpendicular to the dipolar field-lines. 


(A6) 


Appendix B 
Locations for Visible Emission 


The visible pulsar radio pulse emission is modeled in a 
geometry, in which the emission originates from dipolar field- 
lines within the open-field region, and visible emission is 
emitted tangentially to the field-lines at the source points 
(Cordes 1978; Hibschman & Arons 2001; Kijak & Gil 2003) 
and parallel to the line-of-sight direction. In this geometry, the 
location for the point of visible emission can be expressed in 
terms of the polar and azimuthal angles (pv, pv) in the 
magnetic frame (subscript b) at a particular pulsar phase, ~. For 
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given ¢ and a, the visible point is defined by (Gangadhara 2004; 
Yuen & Melrose 2014) 


cos 20py = =(©0s On 8 + cos? Om — sin? On), 


sin¢ sinw 


, (B1) 
sina cos C — cosa sin Å cos w 


tand,y = 


where cos Ôm = cos a cos Ç + sina sin Ç cos(¢ — Y). By using 
Equation (C5), the (@,y, @p,y) pairs can be transformed to (Ay, 
oy) in the observer’s frame. The values of 6,y and ¢,y change 
as a function of 7, and hence the visible point moves as the 
pulsar rotates and traces a closed path after one pulsar rotation, 
referred to here as the trajectory of the visible point. For our 
simulation described in this paper, the relevant solutions are 
obtained from the nearer of the two magnetic poles to the line-of- 
sight, relative to Y = 0°, and the impact parameter, 8 = ¢ — a, is 
minimum. 

It is difficult to determine the height of a radio emission 
source, but it is usually estimated from (i) relativistic phase 
shift or (ii) geometry. The former model requires an 
asymmetric pulse profile that possesses a clearly discernible 
cone and core components (Dyks et al. 2004), whereas the 
latter model specifies the emitting source being located on the 
last closed field-lines (Kijak & Gil 2003). In this paper, we 
adopt model (ii). The assumption that visible emission only 
originates from within the open-field region suggests the 
dependence of the source points on height. There exists a 
minimum height, ry for each w where the trajectory, for a given 
Ç and a, is tangent to the locus of the last closed field-lines, 
defined by Yuen & Melrose (2014) 


— FE sin? Oby 
sin? Os (b,y) sin Or (dy) 


(B2) 


Here, the angle between the point where the last closed field- 
line is tangent to the light cylinder and the rotation axis is 
designated by 6,, and 0p; represents the polar angle of the point 
on the last closed one. In addition, ry reaches its minimum and 
maximum at y = 0 and w = 180°, respectively. Equations (B1) 
and (B2) imply that visible emission is from only one point in 
the magnetosphere for any given y, denoted by (ry, 0v, dy). In 
this model, the center of the profile window is located at 
w= 0°, where the line-of-sight, and the magnetic and rotation 
axes are coplanar. This implies that emission is from lower 
(higher) altitude at the center (edges) of a profile (Gangadhara 
& Gupta 2001; Karastergiou & Johnston 2007). Emission from 
highly relativistic emitting particles with a large Lorentz factor, 
y, is confined to a narrow forward cone giving a half angle of 
size 1/y to the magnetic field line. The angle is small due to 
aberration resulting in some spread in emission about the 
direction tangent to the field-line. Therefore, observable 
emission is emanated from around a small range of field-lines 
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inside the last closed field-line, which corresponds to a small 
range of heights r > ry. 


Appendix C 
Coordinate Transformations 


We assume the Cartesian coordinate system, in which the 
arrangement of the rotation and magnetic axes of a pulsar is 
such that Z = &, and Z, = m, respectively. The corresponding 
unit vectors can then be represented by £ĉ, §, Z and X,, §,, Zp. 

The transformation between the unit vectors is established 
through the following 


Xp £ £ Xp 
J, | = R| ĵ | and | ĵ | = RT| ĵ, |, (C1) 
ĉp z 2 2p 
and we have, 
cosacosy cosasiny —sina 
R=] -—sinw cos Y O |, (C2) 
sinacosw sinasinw cosa 


where R" signifies the transpose of R. Similarly, the unit 
vectors for the radial, polar and azimuthal directions in 
spherical coordinates can be described by fF, 6, È and 
P, Ô,, Py» with the transformation being given by 


j ŝ 
ô|=P|$], (C3) 
Q ĉ 
with 
sindcos@ sinsin cosé 
P = | cosĝ coso cos@sing —sinð |. (C4) 


—sing cos @ 0 


The transformation between the angles in the magnetic and in 
the observer’s frames is given by 


cos 0 = cos a cos 0p — sina sin 0, cos ¢,, 


tan(g — Y) = 


sin 6; sin @, (C5) 


cos a sin 0p cos ġ, + sina cos 6, 
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